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1. Introduction 

Understanding whether exotic systems do or do not exist is one of the few 
remaining challenges for the Standard Model. Such states are those with ex- 
plicit gluon, <?, degrees of freedom and/or quarks, q, in color states that are not 
singlets. Hadrons with unconventional quantum numbers, i.e. those not possi- 
ble in qq or qqq systems such as J PC = 1 h , may or may not be exotic. For 
example, one of the key predictions of the Coulomb gauge model described in 
this paper is that the observed 1 h 7Ti(1400) is not an exotic but rather a more 
conventional meson-meson resonance with quarks in two color singlet states. 
Addressing the above challenge has motivated the development and application 
of this model which provides a robust, relativistic field theoretical approach 
that retains the attractive, insightful wavefunction picture. Since QCD is in- 
trinsically a non-perturbative many-body problem, an established, successful 
nuclear structure method has been adopted which entails three elements: 1) 
an effective Hamiltonian; 2) a truncated model space; 3) solving the equations 
of motion using standard many-body techniques. The Coulomb gauge model 
therefore uses several many-body techniques, Bardeen-Cooper-Schrieffer [BCS], 
Tamm-Dancoff [TDA], random phase approximation [RPA], coupled channels, 
variational and exact diagonalization to solve the equations of motion for an 
approximate QCD Hamiltonian in a model space spanned by a truncated num- 
ber of particle Fock states. For this method to work the Fock space expansion, 
|^ >= \qq > +\gg > +\qqg > + \qqqq > +■■■, must involve dressed (constituent) 
quaisparticle states, and not bare (current) quarks and gluons. This is because 
the dressed partons have much larger masses which, from energy considerations, 
enables a sensible truncation to describe the low mass hadron spectrum in anal- 
ogy to the reasonable 1 particle-1 hole, 2 particle-2 hole, etc. truncations for 
low lying nuclear states. 

The following sections detail the Coulomb gauge model beginning with the 
Hamiltonian in section 2 and the BCS vacuum treatment for dressing the bare 



quarks and gluons in section 3. Applications to hadrons are discussed in sec- 
tion 4 for two-body qq (mesons) and gg (glueballs) states, section 5 for three- 
body ggg (oddballs) and qqg (hybrids) states, and section 6 for four-body qqqq 
(tetraquarks) states. Section 6 also addresses dynamic mixing of two and four- 
body states. Finally, a summary is presented in section 7. 

2. QCD Coulomb gauge model 

In the Coulomb gauge the exact QCD Hamiltonian is 

Hqqd = H q + H g + H qg + He (1) 
H q = J dx*t(x)[-ia- V + H*W (2) 

H g = iydx[j- 1 n (x)^ (x) + B (x)-B°(x)] (3) 

H qg = g ydxJ a (x)-A a (x) (4) 

H C = -y J dxdyJ- 1 p a ( X )K ab (x,y)Jp b (y) , (5) 

where g is the QCD coupling constant, \t is the quark field with current quark 
mass m, A a = (A a ,Ag) are the gluon fields satisfying the transverse gauge 
condition, V-A a = (a = 1, 2, ...8), II a = — E° r are the conjugate momenta, 

E? r = -A a + s (l-V- 2 VV)f^A c (6) 
B a = V x A a + X -gj ahc K h x A c , (7) 

are the non-abelian chromodynamic fields, J° = ^ (x)aT"f (x) are the quark 
currents, p a (x) = f t (x)T a *(x) + / a6c A b (x) • IF(x) are the color densities 
with standard SU(3) color matrices, T a — 4p, and structure constants, f abc . 
The gauge manifold curvature is measured by the Faddeev-Popov determinant, 
J = det(.M), of the matrix, Ai = V • D, with covariant derivative, D afc = 
5 ah V-gf abc A c , which governs the kernel, K ab {x,y) = (x, ^M^V 2 M'^y , b). 
In this gauge, the color form of Gauss's law, which is essential for confinement, 
is satisfied exactly and can be used to eliminate the unphysical longitudinal 
gluon fields. The Hamiltonian is renormalizable, preserves rotational invariance, 
avoids spurious retardation corrections, aids identification of dominant, low en- 
ergy potentials, permits resolution of the Gribov problem and has only physical 
degrees of freedom (no ghosts). The standard normal mode field expansions are 



(bare quark spinors u, v, helicity, A = ±1, and color vectors £0=1,2,3) 

* W = /(^3[^(k)6Ac(k) + «A(-k)4 c (-k)] e 4k - x ec (8) 

A " (x) = /(^71I [aa(k) + aat(_k)]eik ' x (9) 

^^-[a a (k)-a a t(- k )K k - x (10) 
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Here &Ac(k), d\c{— k) and a„(k) (/i = 0, ±1) are the respective quark, anti-quark 
and gluon Fock operators, the latter satisfying the transverse commutation re- 
lations 

[a£(k),<#(k')] = (2^) 3 5 ab «5 3 (k-k')[V-(-l) Mfc ^] 

= (2^) 3 <5 ab <5 3 (k-k')^(k) , (11) 

due to the Coulomb gauge condition, k • a a (k) = (— l) M fc M a" (k) = 0. 

The Coulomb gauge model Hamiltonian, Hcg, is obtained by replacing the 
Coulomb kernel with a calculable confining potential and using the lowest order, 
unit value for the the Faddeev-Popov determinant 

H CG = H q + H GG + H G g G + H GG (12) 

Hg G = \j dx I na ( x ) • na ( x ) + B "( x ) ' B °( x )] ( 13 ) 

Hg G = \ jd*dyJ? (x) (j w - ^) C/(|x-y|)j;(y) (14) 

i/CG = -IJ d xdyp a ( X )V(\ X -y\)p a (y) . (15) 

Also, using Maxwell's equation, the quark-gluon Hamiltonian has been con- 
verted to a quark-quark hyperfine interaction involving an effective massive 
(m g = 600 MeV) gluon exchange potential (modified Yukawa kernel) 



ln-°- 62 (^+0.82) 

' i;: "" : p>m g 



U(p) = < P ln°-»(j£+1.41) 



5-509 < 
p 2 +m 2 P ^ mg 



(16) 



Confinement is provided by the Cornell potential, V(r) = —a s /r + ar, or in 
momentum space, V(p) = — 47ra s /p 2 — 8tt<t/p 4 . The string tension, a = 0.135 
GeV 2 , and a s — 0.4 have been independently determined and set the scale for 
the calculation. The remaining parameters are the bare (current) quark masses, 
m u = m d = 5 MeV, m s = 80 MeV, m c = 640 MeV and m b = 3.33 GeV. 
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3. BCS vacuum and gap equations 

The Bardeen-Cooper-Schriffcr method is used to obtain the ground state 
(vacuum). This entails rotating the field operators 



B\c (k) = cosy6 AC (k)-Asiny<4e(- k ) 



Dxc(-k) = cosyd A c(-k) + Asiny&^ c (k) 

c* a (k) = coshe fe a a (k) +sinhe fe a at (-k) , (17) 



producing the dressed, quasi-particle operators a a , B\c and D\c- The BCS 
quasi-particle vacuum, |f2), is defined by B\c\ty — D\c\£l) = = 0, and 

builds on the bare parton vacuum, |0), b\c\0) = d\c\0) = a ^|0) = 0, 

K uark ) = [e~ /(^ AtBn ^ 6 Ic( k )4c(-k) ] )0) (18) 

\n gluon ) = [e"/^^^^'^)^'^^^ |Q) (19) 

Quark and gluon condensates (correlated qq and gg Cooper pairs) naturally 
emerge in the composite BCS vacuum, |f2) = \ft quark ) <S> \^ g iuon)- A variational 
minimization of the model ground state, 8(U\Hcg\^) = 0; yields the constituent 
quark and gluon gap equations 

ks k ~m q c k = | J -^0^ [V(\k - q\)(s k c q x - s q c k ) 

-A(c kSq U(\k - q|) - c q s k W(\k - q|))] (20) 

2 ,2 3 / <iq .,„„ (l + x 2 )(u|-^) 2 (l-i 2 ), 



- * - i y i#i v '<i k - ^ — - • (22 » 

with Sfc = sirufrh, c k — cos<f> k and x = k • q. Here = </>(&) is the quark 
gap angle related to the BCS angle 9 k by, tan(4> k — 8 k ) = m/k, and to k = 
ke~ 2Qk is the effective gluon self energy. The quark gap equation is UV fi- 
nite for the linear potential, Sna/p 4 , but not for the Coulomb potential, 
— 4ira s /p 2 . The gluon gap equation has both logarithmical and quadratical 
UV divergences and an integration cutoff, A = 4 GeV, determined in previ- 
ous studies is used in both equations. The gap equations yield the quark, 
E k = y/M{k) 2 + k 2 = M{k)/sin4> k , and gluon, uj k , self-energies from which 
the dressed (constituent) quasi-particle masses can be extracted at zero mo- 
mentum: M„(0) = 125 MeV, M g (0) = w(0) = 800 McV. The gap angles 
also determine the quark, (qq) = (fi|*(0)*(0)|fi) = -(177MeV) 3 , and gluon, 
{aGp V G£ v } — (433 MeV) 4 , condensates which are in reasonable agreement with 
QCD sum rule, -(236 MeV) 3 , and lattice, (441 MeV) 4 , values, respectively. 
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4- Meson and glueball states 

Turning to excited states, the lightest in the quark and gluon sectors are two- 
body mesons (qq) and glueballs {gg). The meson states were calculated using 
both the TDA and RPA, the latter being essential to preserve chiral symme- 
try. The TDA and RPA states are respectively, \^tda) = <3ljpc( TZM )l fi )' 
\^rpa) — Q n JPc(RPA)\nRpA), with J, the total angular momentum, P, the 
parity, C, the C-parity, n, the radial-node quantum number and 

Ql JPC (TDA) = J ^*% p C (k)^t( k )4t ( _ k) (23) 

Q\ lJPC {RPA) = l^[X^ pc B^0<.)Dl\-k) - Y^ PC B c ^)D%-n- (24) 

Note the RPA now involves the RPA vacuum defined by Q n jPc{RPA)\Q,RpA) — 
0. The TDA and coupled RPA equations are given in Refs. Q3I2]- For the vector 
mesons p, J/tp, the wavefunctions have been generalized to included both s and 
d waves leading to four coupled RPA equations detailed in Ref. [3]. 

Exactly solving the TDA and RPA equations of motion, H CG \^ nJPC ) = 
M n jPc\^ nJPC ), yields predictions for the u/d, s, c and b meson spectra in rea- 
sonable agreement with observation. Several results are especially noteworthy. 
Except for the pion and light j], the TDA and RPA spectra were essentially 
identical. For light mesons, a realistic pion mass (« 150 MeV) was obtained 
along with Regge trajectories that were also consistent with scattering mea- 
surements, J = a(t) = bt + a(0) « M + .5, t = M^ JPC . For heavy mesons, 
the predicted small hyperfine splittings were in excellent agreement with data 
as summarized in Table 1. Further, the CG model accurately predicted [3] the 
lightest bottomonium state years before its discovery [4] and much closer than 
both lattice [5] and non-relativistic QCD [5] results. 



Table 1: Model comparison of heavy meson predictions with data. 



in MeV 


NRQCD 


lattice 


Coulomb gauge 


data 




104 


90 


125 


117.7 




39t 


61+ 


70+ 


71.4 




9421 ± 11^ 


9409+ 


9395 


9389 



predicted before pseudoscalar bottomonium discovery 



It is significant that the CG model is able to simultaneously describe the 
small charmonium/bottomonium and large n/ p splittings with the same hyper- 
fine interaction. This is possible because the RPA meson operator commutes 
with the chiral charge, Q 5 = J dx^ (xjj^ (x) , and preserves chiral symmetry 
yielding a light Goldstone pion. Conversely, the TDA operator does not com- 
mute and the TDA pion for the same Hcg bas mass about 500 MeV. Hence 
chiral symmetry is predominantly responsible for the large ir/p splitting, con- 
tributing about 400 MeV. Note that the model also describes the excited state 
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spin splittings, typically less than 200 MeV, which are entirely from the hy- 
perfrne interaction. These states are not governed by chiral symmetry and the 
TDA and RPA results agree to within a few percent. See Ref. [3] for a more 
detailed discussion of model meson results. 

In the gluonic sector, the TDA glueball wavefunction for two constituent 
gluons is given by 

l^? C > = / j^*i P sZx^)<M«ll(-m) ■ (25) 

Here L is the orbital anglular momentum and S — 0, 1 or 2 is the total in- 
trinsic gluon spin. Using this quasiparticle basis the excited glueball spectrum 
was computed by diagonalizing the Hamiltonian in the TDA truncated at the 
lp-lh quasiparticle level. The RPA spectrum was also calculated and agreed 
with the TDA to within 1 per cent. The predicted J PC = ±+ ,2 ±+ and 3++ 
states agree well with quenched lattice gauge measurements. The predicted 
J PC = 2 ++ ,4 ++ glueballs yielded a Regge trajectories close to the observed 
pomeron result, ap s» .25i + 1. See Ref. [7] for further details. 



5. Oddball and hybrid meson states 

Predictions for the low lying spectra of ggg glueballs and qqg hybrid mesons 
are now presented and discussed. Since these hadrons consist of 3 constituents, 
the masses for selected J PC states are computed variationally 

M -'pc - (yjpc^jpc) ■ W 

The variational approximation has been comprehensively tested in two body 
systems by comparison with exact diagonalization and found to be accurate to 
a few percent. For C = — 1 glueballs (oddballs), Fock states with at least 3 
gluons are necessary and the wavefunction is (<jj=i 2,3 are cm gluon momenta) 

l* JPC > = / d(adq 2 dq 3 S(q 1 +q 2 + ci 3 )F^ fl3 {(i 1 ,(i 2 ,< l3 ) 

C abc <( qi)a H (q2)a ct (q3)|0) (27) 

The color tensor, C abc , is either antisymmetric, f abc , for C = 1 or symmetric, 
d abc , for C = — 1 and Boson statistics requires C = — 1 oddballs to have a 
symmetric space-spin wavefunction. Using a two-parameter variational radial 
wavefunction, the J oddball states have been calculated. The hyperfine in- 
teraction, H qg , was suppressed and only the Abelian component of the magnetic 
fields was included. The Monte Carlo method with the adaptive sampling algo- 
rithm VEGAS was used and numerical convergence required between 10 5 and 
10 6 samples. In Table 2 results [5] are compared to lattice gauge results [SI [TU] 
and a Wilson- loop inspired model [TT]. The oddball mass sensitivity to both 
statistical and variational uncertainties was a few per cent. 
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Table 2: Oddball quantum numbers and masses in MeV. Error (Monte Carlo) 
for H C g is less than 100 MeV, lattice errors are typically 200-300 MeV. 



Model 


1— 3— 5— 7— 


Coulomb gauge [8] 
lattice [9] 
lattice [TD] 
Wilson-loop [H] 


3950 4150 5050 5900 
3850 4130 
3100 4150 
3490 4030 



The predicted oddball Regge trajectories from several approaches are dis- 
played in Fig. 1. Constituent gluon predictions are represented by boxes, solid 
triangles and solid circles and correspond to a Wilson-loop inspired potential 
model [11] , a simpler harmonic oscillator calculation [8] and the Coulomb gauge 
model 8J, respectively. Lattice results are depicted by open circles [9] and dia- 
monds [10 . The odderon trajectories for the harmonic oscillator and Coulomb 
gauge models are represented respectively by the solid lines, ckq = 0.18^ + 0.25 
and q§ g — 0.23t — 0.88, while the uj trajectory is the much steeper dashed line. 



Oddball Chew-Frautschi plot 

C=P= -1 pure glue states 




Morningstar and Peardon 
Meyer and Teper 
Kaidalov and Simonov 
co trajectory 
H M (this work) 
Coulomb gauge H _ (this work) 



40 



M 2 (GeV 2 ) 

Figure 1: Comparison of the u) meson Regge trajectory to odderon trajectories 
from constituent gluon models and lattice. 
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Three key results follow which are also supported by a more recent, but 
simpler constituent gluon model 12 . First, the oddcron starts with the 
state and not the 1 which is on a daughter trajectory. Note that there are 
no lattice 5 glueball predictions which are necessary to confirm this point 
and we strongly recommend that future studies calculate higher J states. 
Second, all approaches predict the 3 mass is near 4 GeV. Third, the predicted 
odderon has slope similar to the pomeron but intercept clearly lower than the u 
value. This provides an explanation for why the odderon has not been observed 
in total cross section data which is predominantly governed by the larger ui 
intercept. Future searches should therefore focus on differential cross section 
measurements, da/dt, at large t where the odderon trajectory dominates. 

Turning to the hybrid meson system, the non-abelian magnetic field terms 
are now included. For a qqg hybrid the color structure is given by SU C (3) 
algebra, (3 <g> 3) <g> 8 = (8 <g> 8) © (8 <g> 1) = 27 © 10 © 10 © 8 © 8 © 8 © 1. Note 
to obtain an over all color singlet the quarks must be in an octet state, like the 
gluon, which produces a repulsive qq interaction, confirmed by lattice at short 
range, that raises the mass of the hybrid meson. Denoting the momenta of the 
dressed quark, anti-quark and gluon by q, q and g, respectively, the hybrid cm 
system wavefunction is 



= dqdqdg5(q + q ■ 



g)$^(q,q,g)T c a c -i3i c (q)4 c -(q)^t(g)|0).(28) 



With two variational parameters, the hybrid mass was computed using the 
Monte Carlo method which required about 50 million samples for convergence 
with error around ± 50 MeV. The predicted [13] low lying mass spectra for light 
hybrid mesons, with both conventional and unconventional (labeled exotica) 
quantum numbers, are presented in Fig. 2. Note the isospin splitting due to 



M(GeV) 
3.0 

2.5 
2.0 



exotica 



uug (1=0) 



exotica 



4o ++ , 2 + 
^1+- 



uug (1=1) 



.0" 

' l 



Figure 2: Low lying isoscalar and isovector uug spectra. 



quark annihilation (only in the I — channel) which increases the hybrid mass. 
The most significant result is that all hybrid masses, especially the lightest exotic 
1 h state, are above 2 GeV. As summarized in Fig. 3, this is consistent with 
lattice [HI [TS] and Flux Tube model \Tjg\ results and strongly suggests that the 
1 h 7r(1600), and more clearly observed 7r(1400), are not hybrid meson states. 

udg Hybrid Meson 



M(GeV) 

2.5 
2.0 



1.5 



PC 

(l)J =(1)1 



l7Ii(1400) 



Bag Flux Lattice Coulomb Experiment 
model tube QCD gauge 
model 

Figure 3: Confrontation of different hybrid meson models with data. 



The strange ssg and charmed ccg hybrid spectra are illustrated in Fig. 4. 
Due to the hyperfine interaction, the ccg spectra has a slightly different level 
ordering than the uug and ssg results. The ssg and ccg exotic 1 ^ states are 
also in reasonable agreement with both lattice and Flux Tube results. 



M(GeV) 
3.0 

2.5 

2.0 



■o 

i 



exotica 



ssg 



■ i 
■o~ 



5.0 



4.5 



4.0 



exotica 



3 
"1 



.0 

T 



ccg 



Figure 4: Selected low lying ssg and ccg spectra. 
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6. Tetraquark states 

As illustrated in Fig. 5, the SU C (3) color algebra for four quarks produces 
81 color states, 3 fg) 3 (g) 3 <gs 3 = 27 ® 10 ® 10 ® 8 ® 8 ® 8 ® 8 ® 1 ® 1. Two 
are color singlets that can be obtained in four different ways depending on the 
intermediate color coupling: singlet scheme, non-exotic meson-meson molecule, 
and three exotic atoms, one octet and two diquark schemes involving the triplet 
and the sextet representations. The two color singlets, 5 ClC2 5 C3C4 and S ClCi 5 C3C2 , 










W 


?7 • "C • 13-8-8-1 


(6\ 9 


1 (molecule) 


© 

Is/® 



Figure 5: Color singlets from four different representations. 

are linearly independent and form a color space as depicted in Fig. 6 with the 
first along the horizontal axis and the second, in the limit of large N c , vertical. 
For physical N c = 3 they are not orthogonal, as the second is now rotated with 
respect to the first held fixed. However they still span the entire color space 
so that any of the four schemes can be represented as a linear combination of 
the two singlets. This means that the color degree of freedom does not forbid a 
tetraquark transition into two mesons. 



Nc=12 Nc=6 

Nc=3 




Figure 6: Color space. 
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Denoting the momenta of the quarks by qi , q3 , and those of the anti-quarks 
by q2, q4, the cm tetraquark wavefunction is 

|* JPC ) = yrfqirfq 2 rfq3dq4^(qi +q2 + q3 + q4)^f A ^A 3 A 4 (qi,q2,q3,q4) 

^c^i 1 c 1 (qi)< C2 (q2)< C3 (q3)<c 4 (q4)l«> , (29) 

where the color elements, -R^cv depend on the specific color scheme chosen. 
Contributions to the Hamiltonian expectation value involve 4 self-energy, 6 scat- 
tering, 4 annihilation and 70 exchange terms each of which can be reduced to 
12 dimensional integrals that are evaluated in momentum space. Again, these 
were computed |17j by performing Monte Carlo calculations (typically 50 mil- 
lion samples) and the hyperfine interaction was not included. The meson-meson 
molecule yields the lightest mass state for a given J PC which is due to cancel- 
lation of certain interactions by color factors in the singlet-singlet molecular 
representation and also the presence of repulsive forces in the other, more ex- 
otic, color schemes. The ground state is a non-exotic 1 ++ with mass around 
1.2 GeV. The remaining low lying spectra for states having both conventional 
and unconventional (loosely labled exotica but not exotic) quantum numbers in 
the molecular singlet color representation is displayed in Fig. 7. As in the hy- 
brid calculation, there are isospin splitting contributions, up to several hundred 
MeV, but only in the octet scheme (not shown) from quark annihilation inter- 
actions (qq — > g — > qq) in the I q q — channel. The annihilation interaction is 
repulsive, yielding octet states with 1 — 2 lower than the I — 1 which are lower 
than the 1 = 0. The meson-meson states are all isospin degenerate producing 
several molecular tetraquark states with the same J PC in the 1 to 2 GeV 



'Singlet color scheme } 

exotica 



-0,1 

- 2 + ~ 
-0 ++ 

- l + ~, ' + 



1=0 1=1 l=2 

Figure 7: Non-exotic tetraquark color singlet (molecule) spectra. 
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region. The I — 1 and 2 states can be observed in different electric charge 
channels (different I z ) at about the same energy, which is a useful experimental 
signature. Most important, the predicted lightest 1 h is at 1.32 GeV, near the 
observed 7r(1400), suggesting this state has a non-exotic meson-meson resonance 
molecular structure. Indeed, the predicted tetraquark mass for all exotic 1 ^ 
states in the octet color configurations is above 2 GeV. This is consistent with 
the model predictions, discussed in the previous section, for exotic hybrid meson 
1 h states also lying above 2 GeV due to repulsive octet color quark interac- 
tions. Note that C parity forbids exotic diquark ( [3 <g> 3] <g> [3 ® 3] ) states in the 
1 h channel. The other J PC states in both the triplet and the sextet diquark 
color representations also have computed masses (not shown) heavier than in 
the singlet representation and comparable to the octet scheme results. Since 
there are two linearly independent color configurations, in lieu of a color mixing 
calculation that would order the eigenvalues for a given basis, the variational 
procedure requires accepting as the best approximation to the physical states 
the molecular scheme as the lightest along with the heavier linear combination 
that is orthogonal. 

Finally, the important issue of meson and tetraquark mixing is addressed |18) 
for the J PC — ±+ and 1 states. Using the notation, \qq > and \qqqq > for 
Y^JPC > ^ -|- ne m j xec i state is given by | J PC ) = a\nn) +b\ss) +c\nnnn) +d\nfiss}, 
where nn = ^{uu + dd). The state \ssss) is not included since its mass is 
much higher. The coefficients a,b,c and d are determined by diagonalizing the 
Hamiltonian matrix. The off-diagonal mixing element, M = (qq\H^ G \qqqq) , 
only involves Hq which connects meson and tetraquark states. Of the six off- 
diagonal matrix elements two, (ss\H GG \nn) and (ss\H GG \nnnn) , vanish and 
one, (nnnn\H GG \nnss) , is computed very small. The remaining three elements 
are, (nfi\H GG \nhnn) , (nh\H GG \nhss) and (ss\H^ G \nnss) . Because of color 
factors, nonzero mixing only exists for qq annihilation between different singlet 
qq clusters. There are two contributions to each mixing matrix element. One is 

Mi = lJ d qi dq 2 dq 3 V(k)Ul ( qi )U K (-q 4 ) (30) 
<(q 3 )VA 2 (q 2 )$ A 7 f5 3A4 (q 1 ,q2 I q3)^fA ( :(-2q 4 ), 
with V(k) the confining potential, q 4 = qi — k, k — q 2 + q3 and 

<31) 

are dressed, BCS spinors. Again s q = sin 4> q is obtained from the gap equation 
solution. The other contribution has the form 

M 2 = IJ dqidq 2 dq 3 U(fc)Vl 4 (q 4 )VA i (-q 1 ) (33) 

U{ 3 (q 3 )V A2 (q 2 )$S 3A4 (qi, q 2 , q 3 )$^(2q 1 ) . 
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Figure 8: Unmixed and mixed f spectrum compared to data. 



Because new model masses are computed, the unmixed variational basis states 
need not be ones producing a minimal, unmixed mass. This allows adjusting one 
of the two variational parameters, denoted by 7, to provide an optimal varia- 
tional mixing prediction. For ++ states, the mixing term vanishes for 7 = and 
then increases with increasing 7. For the optimum variational value, 7 — 0.2, 
the matrix elements are (ss\H£ G \nnss) = 365 MeV, (nn\H^ G \nnnn) = 166 
MeV and (nn\H^ G \nnss) = 45 MeV. With the calculated matrix elements 
and unmixed meson and tetraquark masses, the complete Hamiltonian matrix 
was diagonalized to obtain the expansion coefficients and masses for the corre- 
sponding eigenstates illustrated in Fig. 8. Mixing clearly provides an improved 
description for the /o spectrum as the a meson mass is lowered from 848 MeV to 
776 MeV, the strange scalar meson mass decreases from 1297 MeV to 1006 MeV, 
closer to the observed value of 980 MeV, and the other /o states are now also in 
better agreement with data. Including chiral symmetry corrections, which are 
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omitted by the TDA variational basis states, will further lower the a predic- 
tion, closer to the now accepted value of 450 MeV [T!5] . The mixing coefficients 
also provide new structure insight, predicting the cr//o(600) is predominantly a 
mixture of nn and nfvnn states while the /o(980) is mainly ss and nfiss states. 
This is consistent with the growing consensus that the cr// (600) state is a irn 
resonance (pole in the tttt scattering amplitude) with a molecular tetraquark 
nature. Future mixing calculations will include both scalar glueballs and hybrid 
mesons which should further improve describing the high lying /o spectrum, 
anticipated to have a few newly discovered states, and also aid identification of 
gluonic states. 

For h states, the value 7 = 0.5 yields reasonable n and rf masses with 
mixing elements (nfi\H^ G \nhnn) — 219 MeV, (nh\H GG \nnss) = 157 MeV and 
(ss\H GG \nnss) — 138 MeV. The unmixed 77, rf masses changed from 610 MeV, 
1002 MeV to 531 MeV, 970 MeV, respectively, both closer to the observed values 
of 547.51 MeV and 957.78 MeV. 

A novel mixing result was obtained for the 1 states. Again the mixing 
matrix elements were for 7 = but, and very interesting, also essentially 
for all values of 7. The Coulomb gauge model therefore predicts minimal flavor 
mixing for vector mesons which agrees with the known, predominantly ideal, 
Lu/(p mixing. Related, the model still provides a good vector meson spectrum 
description since the unmixed nn and ss states were already in agreement [3] 
with observation. 

7. Summary 

Concluding, the Coulomb gauge model provides a comprehensive, unified 
quark-gluon framework for realistically describing the vacuum and meson spec- 
trum and also agrees with glueball and hybrid meson predictions from alterna- 
tive approaches. The model is sufficiently robust, as evidenced by accurately 
predicting the rjb mass and attending hyperfine splitting, to guide experimen- 
talists in future particle searches, especially states with explicit gluonic degrees 
of freedom. Further, this approach retains the attractive wavefunction picture, 
not available through lattice QCD, which provides deeper hadronic insight that 
should be helpful in understanding states with unconventional quantum num- 
bers like the observed tti which appears to be a meson-meson resonance. Finally, 
the approach is amendable to further refinements through improved confining 
interactions and extended model spaces involving additional quasiparticle Fock 
states. 
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